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Selected results from r analyses performed using the BABAR detector at the SLAC National Accelerator Lab- 
oratory are presented. A precise measurement of the r mass and the t+t - mass difference is undertaken 
using the hadronic decay mode — > ■k + tt~ tt^Vt . In addition an investigation into the strange decay modes 
r~ — ♦ Kgir~ir°v T and r~ — » A*g7r~ ' v T is also presented, including a fit to the t~ — * Kgn~u T invariant mass 
spectrum. Precise values for M(K*(892)) and T(K*(892)) are obtained. 



1. Introduction 

Although the BaBAR detector [l| was conceived as 
an experiment to test CP-violation in the B meson 
system, the cross-section for r pairs (0.9 nb) is almost 
as high as that for 66(1.1 nb). This makes BABAR an 
excellent place to study r physics. I present prelimi- 
nary results from studies of tau decays to hadronic fi- 
nal states which provide precise measurements on the 
t mass, the t + t~ mass difference (Section 2), and in- 
formation on the branching ratios and mass spectra 

K%tt~v t and 



of the strange hadronic decays r 
r~ — > Kgir~ir°v T (Section 3). 



2. Precise measurement of r mass and 



T 1 T 



mass difference 



A key test of CPT invariance is to measure the dif- 
ference in mass between a particle and its antiparti- 
cle. Using 423 fb _1 of data from the BABAR detector, 
a pseudomass endpoint method was used to measure 
the mass of the r lepton Q. The significant advan- 
tage in using this method is that it allows the mass of 
the r + and t— to be measured independently, which 
allows us to test CPT invariance. 

The current world average of r mass is 1776.84± 
0.17 MeV/c 2 '3], and the mass difference is M T + - 
M T - < 2.8 x 10~ 4 at 90% confidence level. 

The pseudomass endpoint method was first used by 
the ARGUS [3] collaboration, and has since been em- 
ployed by BELLE [|. The premise is first to consider 
reconstructing the mass of the r from the final state 
hadronic products (t~ — > h~v T ): 



M T = y/M 2 h + 2(V5/2 - E* h )(E* h - P* cos 9*), (1) 

where Mh, Eh and Ph are the mass, energy, and 
magnitude of the three-momentum of the hadronic 
system h, respectively, 9* is the angle between the 
hadronic system and the v T \ the * represents quanti- 
ties in the e + e~ center-of-mass frame. The relation- 
ship E* = y/s/2 is used, where s/s is the initial e + e~ 
CM energy (10.58 GeV), and E* is the energy of the r. 



In the above representation the angle 9* is unknown as 
the neutrino escapes undetected; we therefore define 
the pseudomass M p with the condition that 9* = 0. 
This simplifies the above equation: 



M p = 



Ml 



■2(V~s/2-E* h )(El-P*). (2) 



The distribution of M v 

at M n 



± p has a sharp kinematic cutoff 
j. p — M T , although there will be smearing due to 
initial and final state radiation, and limited detector 
resolution. 



The signal mode chosen for this study is 
i — ► 7r tt^tt^i't-, due in part to its high branching 
fraction, (9.03 ± 0.06) xlO~ 2 3], and the ability to 
obtain high signal purity. The data are then fitted to 
an empirical function: 



F(x) = (ps +p4x)tan 



-i (Pi - x) 



1>2 



+ p 5 +p 6 x, (3) 



where the pi are the parameters and x is the pseu- 
domass. A relationship between p\ ,the endpoint pa- 
rameter, and the pseudomass is required. This is de- 
termined from Monte Carlo studies; three different 
samples are generated with different tau mass values 
and these are fitted with the function defined above. 
A straight line is then fitted to the resulting p\ val- 
ues, which provides the relation to M T . The data are 
split into two samples according to the total charge 
of the 3tt hadronic final state, and each sample is an- 
alyzed independently. The results of these fits are 
shown in Figure [T] This yields an M T of 1776.68 ± 
0.12(stat) MeV/c 2 . A number of systematic effects arc 
investigated, including potential uncertainties on en- 
ergy loss measurements for charged particles, and un- 
certainties in the magnetic field. However, the dom- 
inant uncertainty is found to be due to an underes- 
timation of the reconstructed track momenta in the 
detector model. This contributes 0.39 MeV to a total 
systematic uncertainty of 0.41 MeV. 



3. Strange hadronic tau decays 

Strange hadronic tau decays offer a very clean envi- 
ronment for studying the weak current. The branch- 
ing ratios feed directly into a measurement of V us , 
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Figure 1: Combined r + and r~ pseudomass endpoint dis- 
tribution. The points show the data, the curve is the fit 
to the data, and the solid area is the background. The in- 
set is an enlargement of the boxed region around the edge 
position showing the fit quality where pi is most sensitive. 



and fits to the mass spectra can yield resonance pa- 
rameter values which can further our understanding of 
the dynamics of these systems. In this section stud- 
ies are presented of the hadronic mass distributions 
for the decays r~ — > Kgir~v T and r~ — > Kg^n ^ 
(throughout the note, charge conjugate modes are im- 
plied). A fit to the invariant mass spectrum from 
t~ — ► Kg~K~v r is presented along with precise res- 
onance parameter values of the dominant K*(892)~. 
Due to this mass spectrum having a peaking back- 
ground from t~ — » Kgir^TT ^, the hadronic mass 
spectra and branching ratio from this mode were also 
measured and the results used directly to improve our 
Monte Carlo modelling. 

3.1. Analysis of r — ► K^it~tt p t 

To select events of the type e + e~ — > t + t~ with 
one tau-lepton decaying to Kgir~ "7r°zy T , the event is 
first divided into two hemispheres in the center-of- 
momentum system (CMS) using the thrust axis. One 
hemisphere of the event is required to contain only 
one charged track; this is defined as the tag hemi- 
sphere. The other hemisphere is required to have 
three charged tracks; this is called the signal hemi- 
sphere. The tag track and at least one of the signal 
hemisphere tracks are required to originate from the 
interaction point. 

Approximately 35% of r-leptons decay to fully lep- 
tonic final states. Requiring the track in the tag hemi- 
sphere to be identified as an electron or muon while re- 
quiring the signal hemisphere to contain only hadrons 
strongly reduces backgrounds from e + e~ — > qq events. 
Electrons are identified using specialized likelihood se- 
lectors, whereas a neural network is used to identify 
muon tracks. 



Kg candidates are constructed from any two oppo- 
sitely charged tracks with an invariant mass within 
25MeV/c 2 of the K s mass, 497.672 MeV/c 2 Q. Only 
events with exactly one Kg candidate are retained. 
The track from the signal side not originating from the 
K s candidate is required to be identified as a pion and 
originate from the interaction point. Pions are iden- 
tified by dE/dx in the tracking system, the shape of 
the shower in the calorimeter and information from 
the DIRC. All tracks on the signal side are required 
to lie within the geometrical acceptance region of the 
EMC and DIRC to ensure good particle identification. 

In addition, the net charge of the event must be zero 
and the thrust of the event is required to be greater 
than 0.85 to reduce the non-r background. 

Backgrounds from Bhabha events are suppressed 
by requiring the momentum of the tag-side track to 
be less than 4.9GeV/c. Backgrounds from radiative 
Bhabha and /x-pair events with a converted photon 
are suppressed by requiring the modulus of the cosine 
of the decay angle to be less than 0.97. The decay 
angle is defined as the angle between the momentum 
of the 7r + originating from the Kg in the Kg's rest 
frame and the K s momentum vector in the labora- 
tory frame. When this quantity is calculated for e + e~ 
conversion pairs misidentified as pions, its value is con- 
centrated near ±1. From studies of missing transverse 
event energy, backgrounds from two-photon events are 
determined to be negligible. We also require exactly 
one identified 7r° in the event; the trajectory of the 
7T° must be within 90 degrees of the K s tt~ momen- 
tum vector. This ensures that the 7r° is more likely to 
be from the same r as the Ksn. The neutral energy 
not attributed to the Kg or the 7r° must be less then 
100 MeV. This should be very small anyway, but the 
cut is to reject unwanted photons. The energy of the 
7T° in the center-of-mass system must be greater than 
1.2 GeV. This cut is to remove the large background 
contribution in the region below 1.2 GeV. Figure. [2] 
shows the distribution of the 7r° energy. 
The branching fraction B(t~ — » K °7r~ir°is T ) is esti- 
mated by 

B(r- - ffOTr-TrV) = 2 l N N ^\ N ^ , (4) 

TT s ig 

where N TT is the total number of t + t~ pairs in the 
data, iV data is the number of selected events in data, 
iV bkg is the number of background events estimated 
from Monte Carlo, and s . is the corrected signal ef- 

! sig o 

ficiency to include Kg and K^ mesons. 

The hadronic mass distributions for the different 
combinations of final state hadrons are also extracted, 
and used to tune our Monte Carlo. Figure [3] below 
shows the mass distribution for KgTr~ from t~ — * 
K s 'ir~Tr v T , overlayed with the new Monte Carlo gen- 
erated using this analysis. 
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CMS PiO Energy [GeV] 

Figure 2: Distribution of the n° energy. 



Table I 8(7 



K°tt ir°v T ) measured in this analysis. 



Sample B(r~ -> jf°7r~7rV) [%] 

e-tag 0.353 ± 0.008 (stat) ± 0.016 (syst) 
/i-tag 0.329 ± 0.008 (stat) ± 0.016 (syst) 
Combined 0.342 ± 0.006 (stat) ± 0.015 (syst) 



3.2. Fit to 



mass spectrum 



The analysis of the decay r~ — > Kgir~v T is a fit of 
the hadronic mass distribution to a parametric func- 
tion describing the resonant structure. From this we 
obtain precise values for the mass and width of the 
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Figure 3: The data (points) and tuned Monte Carlo (blue) 
prediction for the K%n mass distribution from the r~ — > 
KgTT~ '-k°v t signal mode. 



K*(892) as well as information on other resonances 
present in the spectrum. 

We denote the number of events found in bin i 
(without background subtraction) by n,. The predic- 
tion for the expectation value of rij, z/j = J3[fit], can 
be written 



M 

= ^2 / '''-' / '- 

3=1 



(5) 



where Pi is the expected number of background events, 
\ij is the predicted number of signal events in bin j 
before detector effects (the "true" distribution), and 
Rij is a response matrix that reflects the limited ef- 
ficiency and resolution of the detector. The value of 
Rij is 

Rij = P(found in bin i |true value in bin j) , (6) 

and thus the efficiency for bin j is found by summing 
over all bins where the event could be found. The 
predicted number of events in bin j of the true distri- 
bution can be written 



Mj = AHot / /(to; 0) dm , 

Ibinj 



where m denotes the Kgir invariant mass and 



resents a set of parameters. 

The probability density function (pdf) /(to; 



be written 



f(m;0) oc -ll 



(7) 



rep- 



can 



1 



x P \P 2 \F V \ 2 



1 



3{m K 2 -m„ 2 ) 2 2 



where s = 
by 

F v 



1 



1 + /3 + 7 



4s(l + 2^) 
Here the vector form factor Fy is given 

[BW K i(s)+pBW K 2(s)+jBW K3 (s)} . 



This form allows for the K*(892) and two additional 
vector resonances. The quantities (3 and 7 are com- 
plex interference terms between the resonances, and 
the B W terms refer to the to relativistic Breit-Wigner 
functions for the specific resonance, given by 



BW R (s) 



Ml 



s- m 2 r + i^r R ( s ) 

The energy dependent width is given by 

2£+l 



r R ( s ) = r i? 



M% ( P(s) 



where 



a \P(Ml), 
P(s) = ^{ S -M%){s-M 2 _) , 



(8) 



(9) 



(10) 
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and where M_ = M K - M v , M + = M K + M v , and I 
is orbital angular momentum. Thus one has I = 1 if 
the Kit system is from a P-wave (vector), or £ = if 
the Kit system is from an S-wave (scalar). 

The scalar form factor requires a different paramet- 
ric function and can include contributions from the 
K5(800) and K5(1430) signals. This is 



F s = x 



M 2 

if* (800) 



A' 



■BW K * {S00) {s) 
BW K5(im (s) . (11) 



' M 2 

1V1 K*(1430) 

Each of the background modes is subtracted from the 
data, and then a least-squares fit is performed to the 
resulting mass spectrum. The fit model includes Tj 
as a scale factor that relates the luminosity of the 
Monte Carlo sample for mode j to that of the data, 
and factor that allows for the uncertainty in 

the prediction of the rate of the background process. 
The best estimate of rj is equal to unity, but this 
is treated as a Gaussian distributed quantity with a 
standard deviation equal to the relative uncertainty 
on the production rate for the jth background mode. 

The uncertainties in the values of other nominally 
fixed model parameters, e.g., the resonance parame- 
ters of the K*(1410), can be incorporated into the fit 
in a similar way. For a given parameter rj one has a 
previously estimated value f] and standard deviation 
a n , taken from the PDG. One includes in the minimi- 
sation function a Gaussian term in rj centered about 
f] with a standard deviation cr^, and regards rj as an 
adjustable parameter. 

We also include terms in the minimisation which ac- 
count for the uncertainty in the shapes of background 
mass distributions. This is particularly true for the 
KgTr~Kl mode, as it makes a larger contribution and 
the information on its shape is based largely on lower- 
statistics measurements from LEP [l(| • We introduce 
two additional adjustable parameters, a — (01,02), 
which have the effect of shifting and stretching the 
shape of the distribution 11]. This transformation 



is applied to the rriij values for the KqTt 



Kl 



back- 



ground mode and the altered values are then used in 
the minimisation. The fitting procedures described 
above have been carried out using a variety of hy- 
potheses. Figure 4(a) shows that a single K*(892) is 



clearly not enough to model the mass spectrum accu- 
rately. This was seen by the Belle collaboration 
which proposed that the distribution should contain 
contributions from a Kq(800) scalar and K*(1410) vec- 
tor resonances. 



In the region around 1.4 GeV in Fig. 4(a) the data 



are significantly higher than the fitted curve. The addi- 
tion of the K*(1410) gives a significant improvement 
to the high mass region, yielding a x 2 of 130.04 for 
95 degrees of freedom. In these fits the rate of the 
K*(1410) was allowed to vary within the error given 
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Figure 4: A fit to the r — + K%-k v t mass distribution 
using (a) single K*(892) resonance, (b) a combination of 
KS(800) + K*(892) + K*(1410). 



in the PDG. 

The inclusion of the Kq (800) further reduces our x 2 
to 113.05 for 94 degrees of freedom. This is a signif- 
icantly better goodness-of-fit value than our K*(892) 
+ K*(1410) fit model. For the mass and width of the 
Kg (800) we use the measurements from the BES col- 
laboration M — 841 ±301. 73, 
The result is shown in Fig. |4(b)| 

If instead of a K*(1410) one uses a scalar Kg(1430), 
one finds a comparable x 2 value of 114.11 for 94 de- 



618±90l^ 4 . 
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grees of freedom. As such the K*(1410) and Kq(1430) 
cannot be differentiated on their value. To study 
what combination of K*(141Q) and Kq(1430) are 
present, one could exploit the different spins of the 
two resonances by carrying out an angular analysis. 
This is not part of the current study. The resulting 
values for the mass and width of the K*(892) are found 
to be 

M(K*(2Q2)~) = 894.57 ± 049 (stat.) MeV/c 2 
T(K*(892)-) = 45.89 ±0.43 (stat.) MeV/c 

The statistical errors quoted already cover a number 
of systematic uncertainties such as those in the rates 
and shapes of backgrounds, which were incorporated 
by including corresponding adjustable parameters in 
the fit. Several additional sources of systematic uncer- 
tainty are also taken into consideration. The response 
matrix Rij is derived from the Monte Carlo simulation 
of the detector. As a conservative estimate of the un- 
certainty of the detector response, which is dominated 
by modelling of the tracking and Calorimeter, we have 
varied the parameters of the response matrix by up to 
±10%. As a check of the fitting method we have taken 
a fully reconstructed Monte Carlo sample of signal 
events, and fitted them using the signal model. As the 
MC generator models the r~ — > K^tt^v t decay with 
only the K*(892) resonance, the fit model also only 
contained this resonance. A further uncertainty in the 
fit model stems from the choice of resonances. We take 
the difference in the mass and width values for the 
K*(892) between our nominal fit and the alternative 
models which also yield comparable x 2 values, as an 
estimate of the systematic uncertainty. The quadratic 
sum of all these sources of systematic uncertainty lead 
to an error on M(K*(892)) and T(K*(892)) of 0.19 
MeV/c 2 and 0.57 MeV/c respectively. 



K°ir tt u T) which is found to be: 

(0.342 ± 0.006 (stat.) ± 0.015 (sys.))% . 

For the t~ — » Kgir~ir°v T mode we have measured 
the mass distributions of different combinations of fi- 
nal state hadrons: 7r~7r°, Kgir~, KgK and Kgir~ir°. 
These were used to make important improvements to 
the TAU0LA Monte Carlo generator, which allowed for 
a precise estimation of the background contribution 
from this mode in the analysis of the r _ — > Kgir~v T 
channel. 

We have carried out a fit of the hadronic mass dis- 
tribution for t~ — > Kgir~v T . This yields precise mea- 
surements for the mass and width of the K*(892) res- 
onance: 

M(K* (892)~) = 
894.30 ± 0.19 (stat.) ± 0.19 (syst.) MeV/c 2 , 

T(K*(892)~) = 
45.56 ± 0.43 (stat.) ± 0.57 (syst.) MeV/c . 

These values confirm the Belle collaboration's mea- 
surements |8| that indicated a K*(892) mass several 
MeV higher and a width several MeV lower than the 
world average. The results reported here represent 
a factor of two improvement in precision relative to 
the Belle measurements. We analyse the possibility 
of other resonances being present in this mass spec- 
trum, and conclude that a combination of K*(800), 
K*(892) and K*(1410) provides a good description of 
the data. Figure [5] shows the results of various mea- 
surements that went into calculating the 2008 PDG 
average values for the mass and width of the K*(892). 
The Belle 2007 result and our result both indicate a 
shift towards 895 MeV for the mass value. 



4. Summary and Conclusion 

Measurements of the r mass and t + — t _ mass differ- 
ence have been carried out yielding results of: 

M T = 1776.68 ±0.12(stat.)±0.41(syst.)MeV/c 2 . 
MT ~(M) lT+ = ~ 3 ' 4± L3 ( stat ± 0-3(syst.) x 10~ 4 . 

where (M) is the average of M T + and M T -. The r 
mass result is in good agreement with the world av- 
erage. We also find the mass difference result to be 
consistent with the results published by the Belle Col- 
laboration at 1.8(7. 

We have also carried out studies of the decays t~ — > 
K%tt-v t and t~ -> K%tx-tPv t using 384.6fb _1 of 
data. We have measured the branching ratio for r _ — > 
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